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ABSTRACT 


Solar collector is a core sub-system in solar thermal energy systems with various applications including solar drying of agricultural 
products. This research paper studied the effect of input parameters that affect the collector size and efficiency of a solar cabinet 
dryer. For solar drying of tomatoes slices, design equations of collector top heat loss, bottom and edges heat losses, mass of 
moisture removed from crop, quantity of heat required for evaporating moisture from crop, drying chamber heat lost, estimation of 
average solar radiation for location using Angstron-Prescott equation, collector area and collector efficiency were coded in MATLAB 
programming language. The month of August was used as the design month bearing in mind that the month has the least solar 
radiation in Bauchi. Variables studied included single and double glazing, wind speed in the range of 1 to 5m/s, thermal conductivity 
of insulating materials between 0.01 to 0.1W/mkK, collector tilt angle 0 to 60° and initial moisture content of crops in the range of 20 
to 96%. The solar collector area and efficiency were calculated from the MATLAB programme codes. Result shows that wind speed of 
1 and 2m/s have smaller collector areas and relatively higher collector efficiencies compared to wind speeds of 3, 4 and 5m/s for 
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both single and double glazing. Crops with low initial moisture content required small collector areas compared to those with high 
initial moisture content and optimum collector tilt angle of 13.01° for Bauchi. The use of insulating materials of low thermal 


conductivity minimizes heat lost and hence lower collector areas and higher collector efficiency. Collector air gap height of 50mm 
and optimize collector tilt angle of 13.01° in addition to materials of low thermal conductivity for thermal lagging maximizes solar 
energy input utilization, high thermal efficiency and small solar collector area. This design approach minimizes the use of materials, 
will maximized energy usage and give optimum returns on investment. 


Keywords: Collector Area, Optimum Returns on Investment, Parametric Study, Solar Dryer Collector. 


1. INTRODUCTION 


Reduction of post harvest losses significantly contribute to the availability of food. A significant percentage of these losses are 
related to improper and/or untimely drying of foodstuffs such as cereals, grains, pulses, tubers, meat, fish, vegetables etc. Drying is a 
simple process of moisture removal from an agricultural product in order to attain the desired moisture content and can be 
achieved through several methods such as open-sun drying, solar drying and mechanized drying. Solar drying is an energy intensive 
operation and depends on climatic conditions (humidity, wind speed, solar radiation, cloud cover etc) and drying materials’ 
properties (chemical composition, physical structure, size, shape etc) [10]. The removal of moisture prevents the growth and 
reproduction of microorganisms that causes decay, reduces the product weight and volume as well as the transportation cost. The 
drying of agricultural produce under the sun is a common practice most especially in developing countries. Traditional method of 
drying involves spreading the crops on concrete floor, mats, tarred surface roads and other forms of material and turned regularly 
until they are sufficiently dried to the desire moisture content. This method is characterized by a number of shortcomings including; 
lack of process control, non-uniformity in the dried product, soaking by rain, theft and vandalism, contamination by dust, rodents, 
and other domestic animals [9]. Also, prolonged open sun drying often causes deterioration of vital ingredients like vitamins, 
minerals and sensory characters of the dried product and thus, less market value (quality falls below local and international market 
standards). However, this method is more affordable to the farmers. Solar drying method converts sun's radiation into heat for 
drying of agricultural produce at a superior drying rate compared to the open sun drying method with improved quality of dried 
products at a relatively low cost. The solar collector supplies the product with hot airflow either generated naturally due to density 
changes resulting from temperature differences or through forced convection by incorporating a fan or blower in the drying system. 
Solar dryers can be constructed from locally available materials and are useful in areas where fuel or electricity is expensive, land for 
sun drying is insufficient or expensive, sunshine is plentiful but the air humidity is high. Besides, solar drying offers reduction to 
environmental risk compared to the mechanized drying method. Different types of solar dryers with varying sizes and designs have 
been developed and tested with some degree of efficiency. This work investigate the effects and sensitivity of parameters that 
affects collector size and efficiency and, provide a MATLAB programme code that could be used to design for minimum collector 


area required and maximum efficiency for drying. 


2. METHODOLOGY 

Parametric Design of Solar Dryer 

Solar air heater is an integral part of the solar drying system and is basically made up of the glazing, thermal insulation, absorber 
plate, collector air inlet diameter for air flow into the channel and heat storage bed (optional) for improved performance. The above 
units including variables like collector tilt, depth of air passage in collector, wind speed, initial moisture content of crop and mass of 
crop to be dried are vital in determining the size and efficiency of the air collector. It is imperative to study the effect of these 
parameters in determining the size and efficiency of solar collector for improved performance. 


Thermal Insulation 
In solar collectors, the solar energy absorbed by the absorber plate is converted to useful energy gain while part of it is lost as 
thermal losses through the top, bottom and edges of the solar collector. With the assumption that all the losses are based on a 


common mean plate temperature ie , the overall heat loss from the solar collector can be represented as given by [6]: 


Oross a U,A, (T, ~ i ) ooo (1) 
Where; 
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U; is solar collector overall heat loss coefficient, W/m°K; T, is the ambient temperature, K; 
A. is the collector area, m2. 


The overall heat loss is the sum of the top losses, back losses and the edge losses given as: Qhoss = Q, + Q, + Q, ai(2) 


Where; 
subscripts t, b, and e represents the top, back and edge contributions of the heat losses respectively. 


Collector Top Heat Loss for Single Glazing 
The approximate overall top heat loss for the single glazing as proposed by [11]. It assumes that if the mean plate temperature ia 
the ambient temperature To and convective heat transfer coefficient at the glass top h,,, are all known, then the overall heat 


transfer coefficient at the glass top cover U., which includes both the convection and the radiation effects is then computed using 


equation 3. 
-l 
U, = (hp, + Nog ; + (Pepa +h,)" +e 8) 
g& 
Where; 
olf,’ +17,’ \+(r, +7,) 
“= 4 (4) 
== 
é, é, 
12.75[(r, -T, )cos sf? 
cpg (r 1T pes Lo .. (5) 
Pp & PS& 
hp, = 08, (1, +7, \T, +T,) ot 


Where; 


Reng is the radiation heat transfer coefficient between plate and glass, W/m?k 


Ise is the convection heat transfer coefficient between plate and glass, W/m?k 


Ris is the radiation heat transfer coefficient between glass and the ambient W/m?k 


o is Stefan—Boltzmann constant; ie is glass cover temperature, °C 


Li, is air gap spacing between plate and glass, mm 


B is collector tilt angle, ° 


é, is the emmitance of glass cover 

L, is the thickness of glass, m 

k, is the thermal conductivity of glass, W/mk. 

Assuming the sky temperature to be same as that of the ambient temperature, the glass cover temperature, T,, can be computed 


as given in equation 7 by [4]. However, for cases when the sky temperature is different from that of the ambient, the approximation 
proposed by [11] can be used. 
pi l*hth 
i 
oe i 
Where; 


w(f) 
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- li2x10*(r, +0.27,) +h, | +0.31, 
lox10(c, +0.028\r, +0.57,) +0.61,,°°(r, -7, cosa |" 


Ed . (8) 


é, is emmitance of plate 0.1< &; <0.95; ‘i is ambient temperature, °C 


Lie air space between plate age glass cover, m 


h,, is wind heat transfer coefficient (W/m? °C) 0< V <10m/s, if the wind velocity is in V (m/s), then: h,, = 5.7 + 3.8V... (9) 


Collector Top Heat Loss for Double Glazing 
The top heat loss from a double glazing flat plate solar collector can be estimated based on the method also proposed by [11] as 
expressed in equation 10. 


l l 
U, = (per + is + (eres + A vigs ia + (here ot h, a + ane + = wv (10) 
Ko K,> 
Where; 
2 2 
-_ of, Ty : C27, “ay 
—+—-] 
é, Eo} 
12.75|(7, -T,, Cosh? 
cpg) — (r 1T lure 0.21 .. (12) 
P gl pgl 
2 2 
Nygig2 = ot a2 I fe Ta) sil) 
—4f— 1 


é é 


gl g2 

12.75|(7,, -T,. )CospP? 
cpg2 (r LT a 0.21 w+ (14) 
gl g2 


glg2 


o8,.(T..° ry ia Fae +T,) (15) 


h,. 2a — 


Where; 


Poet is the radiation heat transfer coefficient between plate and inner glass, W/m?k 


Noe is the convection heat transfer coefficient between plate and inner glass, W/m?k 


h 


rglg2 |S the radiation heat transfer coefficient between inner and outer glass W/m?k 


Mes is the convection heat transfer coefficient between plate and outer glass, W/m?k 


h 


rg2a_'S the radiation heat transfer coefficient between glass and ambient, W/m?k 


3 and T 5 are the inner and the outer glass temperatures respectively, °C 


L 1 iS air gap space between the absorber plate and the inner glass, m 


Dp 
Lig? is air gap space between the inner and outer glass, m. 


The temperatures of the inner and the outer glasses T and T 9, 


can be approximated according to the method by [4] as 
expressed in equations 16 and 17. 


te ‘ 
on ye ...(16) 
1 
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p Pkt tT 7 
a ey ne 
2 
-1 
Lee NT aT l+e,\T, -T osp > 
sau(t+| M, a) (| vk Te ee A +0.3L,, 
Where; f, = ” =) (18) 
3 0.25 
3A55¢ |T = (2-«, \r, -T,») 0.8 (2-«, \r, -T,, CosB 
Pp Pp 0.2 
6 Lye 6 


lf2x10*(r, +0.27,) +}'+0.32,, (0.7-0.2e, ) (19) 


: lox10-*(c, +0.028\r, +0.57,) +0.6L,., °° (r, -7, Coss}? |' 
Where; 


é, is the emmitance of plate 


Sr 


Bottom and Edges Heat Losses of the Solar Collector 
The thermal losses coefficient for the bottom of a flat plate solar collector as [6]: 


U,=— (20) 


Where; 
kK; is thermal conductivity of the insulating material, W/mk 
X, is thickness of the insulating material, m. 


The edge loss coefficient depends on the area of the solar collector and is usually very small when compared to the losses of the top 
and bottom of the solar collector. The total heat loss coefficient can therefore be given as in the following expression by [6]: 


U,=€,40,70, (21), 
For uniform thickness of the thermal insulation for both edge and back of the collector, then equation 28 becomes U, =U, ): 


U, =U, +2U, (22) 


Amount of Moisture to be Removed from the Product, kg (771, ) 

The moisture to be removed from the product to be dried is given by equation 23 as reported by [11]. 
m,—M, 

nh =n. ... (23) 
100-m, 

Where; 

m_, is the dryer capacity per batch, kg 


M, is the initial moisture content of tomatoes, % 


M ¢ is the maximum desired final moisture content of tomatoes, %. 


Quantity of heat energy required for evaporating moisture from product, Q 
The quantity of heat required to evaporate moisture from the product to be dried is given by [11] as: 


Q=m,, xhj, ... (24) 
Where; 


Q is the amount of energy required for the drying process 
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M1,, is mass of water vapour to be evaporated, kg 


Nye is latent heat of vaporization of water, kJ/kg 
hj, = 4186 (597 —0.56r,, ) 1. (25) 
Where; 


tpr is maximum allowable temperature of the dried product. 


Area of the solar collector unit 
According to [7], the model equation for calculating the area of a solar collector is given as: 


1. Q 
° P,l(ar)I, —U, xtx(Z, -T,)| 


Where; (ca) is the transmittance-absorptance product that represents the effective absorptance of the cover-plate system 


... (26) 


di is total solar radiation incident on the collector, W/m? 


U, is air heater overall heat loss coefficient, W/m2K 
tis the duration, s 


T, is the outlet collector temperature, K 


ii is the collector inlet temperature which is same as the ambient temperature, K 
Q is amount of energy required for the drying process. 


But in strict terms, Q should include the conductive heat lost from the drying chamber, Q. 


Drying Chamber Heat Lost 
The conductive heat lost in the drying chamber is estimated using equation 27 


Q. =U,xtx(T, -T,) ... (27) 


Where; és is the back heat loss from the drying chamber 


Total System Load 


The system total load is given as in equation 28 


OQ, =2+Q. .. (28) 
Solar Collector Efficiency 


The collector's efficiency, 77. is expressed as given by [9] as follows: 


_ pV ,c, AT 
LA, 


1). .. (29) 


Where; 
C , is the specific heat capacity of air, J/kgK 
AT isthe temperature elevation, K 


I. is the solar insolation, W/m? 
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V, is the volumetric flow rate, m?/s 


pis the density of air, kg/m? 


Estimation of Average Solar Radiation 
The most convenient and widely used correlation for predicting solar radiation was developed by Angstrom and later modified by 
Prescott. The Angstrom formula is given by [12]. 


H 
— =at+b— ...(30) 
A 


HZ is the monthly average of the daily global radiation 


fl. is the monthly average of the daily extraterrestrial radiation on a horizontal surface at the location 


n is the monthly average of the sunshine hours per day 


N is the monthly average of the maximum possible sunshine hours per day at the location 
a and b are empirical constants for the location. 


The monthly mean daily extraterrestrial irradiation and monthly mean day length can be derived from the equations 38 and 39 
respectively [4]. 


aoe: 1+0.033| 360—— | |x sin gsin 6| cm, +cos@cosdé sin a, |... (31) 
1 365 360 


oO 


a 2 7 
N =—cos'(- tang tan6) a (82) 
15 
Where; 
the hour angle @,, latitude @ and declination O are related by: 


cos@, =—tan gtand ... (33) 
Where; 


O is the solar angle of declination and is approximately given as: 


5 = 23.45sin 308 (9844 n) .. 34), 
365 


The value 1367W/m? has been recommended for the solar constant ([1]. 
Where; 
n is the day of the year from January 1 to December 31. 


Solution Scheme of the Solar Dryer Design 

The study of the solar dryer parameters for sizing and performance was carried out using MATLAB programming language. Design 

equations for collector top heat loss for single and double glazing, bottom and edges heat losses, amount of moisture to be 

removed from crop, quantity of heat required for evaporating moisture from crop, collector area, drying chamber heat 
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lost, collector efficiency and estimation of average solar radiation were coded in MATLAB programming language. Input data of 


number of glass cover, wind speed, solar radiation available for Bauchi town in the design month of August, thermal conductivity of 


lagging materials, collector tilt angle, collector air depth passage, product design and ambient temperature, crop initial moisture 


content, latitude location and other input data’s were used for the study. To facilitate the parametric design study of the solar dryer, 


refer to Table 1 for flat plate air heating system. 


Table 1 Flat Plate Air Heater System Design Parameters 


Parameter Design value 
Crop Tomatoes 
Mass of crop, kg 6 

Initial moisture content, % 94 

Final moisture content, % 45 
Latitude of Bauchi, degree 7.7306 
Thickness of glazing, mm 4 
Transmissivity of glazing. 0.88 
Thermal conductivity of glazing, W/mk 1.05 

Air gap space between plate and inner glass, mm 10 

Air gap space between inner and outer glasses, mm 10 
Absorber plate (Aluminum) 

Emittance 0.09 
Absorptance 0.91 
Thickness of absorber plate, m 0.05 
Thermal conductivity of plate, W/mk 205 
Specific heat capacity of absorber plate, J/kgK 900 

Back insulation thickness, mm 50 


3. RESULTS AND FINDINGS 
Glazing Number and Wind Speed 


60 


N Ww ~ uw 
Oo oO oO oO 


Collector Efficiency (%) 


ay 
fo) 


Wind speed (m/s) 


Figure 1 Effect of glazing and wind speed on collector area and efficiency 


aq 
£ 
o e=te= N=1 
= N= 
<q 
5 =e N=2 
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S 
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Fig. 1 shows the effect of glazing number and wind speed on collector area and efficiency. Double glazing minimizes collector top 


heat loss and thus, decreased collector area required compared to single glazing signifying more heat retention and utilization 
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(efficiency). Comparatively, wind speed of 1 and 2m/s have low collector areas as against 3, 4 and 5m/s especially for single glazing 


indicating that heat is carried away more faster at high speeds (air flow rate). Thus, more collector area is required to complement 
for the lost heat in the collector for single glazing which transcends to relatively low collector efficiency. For average wind speed of 
2.5m/s obtainable for Bauchi, collector area of 1.88 and 1.63m? are needed for single and double glazing. No significant difference in 
the collector efficiencies for wind speed of 3, 4 and 5m/s for double glazing. 


Glazing Number and Thermal Conductivity of Insulating Materials 

Substantial amount of heat could be loss through the back and edge of absorber plate if not properly or poor insulator is used. 
Lagging of collector systems with materials of synthetic origin of low thermal conductivity are not only expensive (rock wool, fibre 
glass, polystyrene and the polyurethane) but may pose health risk to those handling them in fibrous form. The thermal conductivity 
of insulating materials is proportional to collector area and inversely proportional to the collector efficiency for both single and 
double glazing (see Fig. 2). For this research, collector areas of 1.88 and 1.62m? are required for single and double glazing while, 
efficiency of 45.95 and 53.63% are obtained respectively using saw-dust which is readily available and relatively cheap with thermal 
conductivity of 0.08W/mkK. For any application, the economy viability of using such an insulator is needed. 


70 
60 
5 oe 
~— N 
> ae iS 
§ ; 
‘a 40 o =—@=N=1 
+ <q 
+ 30 o == N=2 
2 g 
o 2 =te=N=1 
> 20 8 
1S) —=<<_=N=2 
10 
0 
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 
Thermal Conductivity of Insulating Materials (W/mk) 
Figure 2 Effect of thermal conductivity of insualting material on collector aea and efficiency 
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Figure 3 Effect of collector tilt on collector area and efficiency 


One of the parameters that affect the performance of a flat plate solar collector is its angle of tilt with the horizontal as the variation 
of tilt angle affects the positioning of the collector to receive maximum amount of solar radiation. As the collector tilt increases from 
0°, the collector area decreases and reaches a minimum at 1.88 and 1.63mé? for single and double glazing corresponding to optimum 
tilt angle of 13.01° and system efficiency of 45.96% and 53.64% respectively (see figure 3). It is worth to note that solar collector 
using monthly average optimum tilt angle will received more maximum solar radiation than that using optimum annual tilt angle. 


Initial moisture content of crop 

Most commonly dried fruits and vegetables have initial moisture content that ranges from 72 to 94% which need to be dried to 
equilibrium moisture content of 7 to 4.5%, while cereals like maize, sorghum have about 30% initial moisture content and needs to be 
dried to a safe storage moisture content of about 14%. From Figure 4, it can be seen that the required collector area is proportional 
to the percentage of initial moisture content of the crop for both single and double glazed collectors. However, double glazed 
collectors require much smaller areas due to the reduction in the top heat loss coefficient as against single glazed collectors. At 94% 
initial moisture content for tomatoes, the collector areas required were found to be 1.88 and 1.62m? for single and double glazing, 
and efficiency of 45.96% and 53.64% respectively (see Fig. 4). 


2.5 


S 
a 


=—@—N=1 


1.0 
ae aa 
os + 


25 30 35 40 45 50 55 60 65 70 75 80 85 90 100 
Initial Moisture Content (%) 


Collector Area (m2) 
= 
un 
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Figure 4 Effect of initial moisture content of crop on collector area and efficieny 


4. CONCLUSION 


For development of an efficient and effective solar drying system, studying some vital parameters that affect the size and efficiency 
of the integral part of the drying system (solar air heater) and the auxiliary heating system for continuation drying operation (if any) 
are essential. The MATLAB programme language could serve as a design format to be use for all crops and locations provided the 
changing variables are known and substituted into the programme. 
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